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EUGEN J. IONASCU 

Abstract. We study the existence of equilateral triangles of given side lengths and with integer 
coordinates in dimension three. We show that such a triangle exists if and only if their side lengths 
are of the form \j2(m 2 — mn + n 2 ) for some integers m, n. We also show a similar characterization 
for the sides of a regular tetrahedron in Z 3 : such a tetrahedron exists if and only if the sides are of 
the form ky/2, for some k € N. The classification of all the equilateral triangles in 1? contained in 
a given plane is studied and the beginning analysis is presented. A more general parametrization 
is proven under a special assumption. Some related questions are left in the end. 



1. INTRODUCTION 

It is known that there is no equilateral triangle whose vertices have integer coordinates in the 
plane. One can easily see this by calculating the area of such a triangle of side length I using the 
formula Area = ^-j^- and by using Pick's theorem for the area of a polygon with vertices of integer 
coordinates: Area = ^ + %i — 1 where %b is the number of points of integer coordinates on the 
boundary of the polygon and %i is the number of such points in the interior of the polygon. Since 
Pick's theorem implies that this area is a rational number of square units, the formula Area = 
says that this is a rational multiple of v3 since I 2 must be a positive integer by the Pythagorean 
theorem. This contradiction implies that no such triangle exists. 

The analog of this fact in three dimensions is not true since one can form a regular tetrahedron 
by taking as vertices the points O(0, 0, 0), A(l, 1,0), B(l, 0, 1) and C(0, 1, 1). It turns out that the 
sides of such regular tetrahedra have to be of the form ky2, k £ N. Moreover and as a curiosity, 
one can use the facts derived in this note to show that there are only three regular tetrahedrons in 
Z 3 having the origin as one of their vertices and of side lengths 9v2, where the counting has been 
done up to symmetries of the cube that leave the origin fixed. In the figure below, that we generated 
with Maple, we show three regular tetrahedrons that together with all their cube symmetries fill 
out the class just described. 



Date: February 2, 2008. 

Key words and phrases, diophantine equations, integers, Quadratic Reciprocity Law. 



1 



2 



EUGEN J. IONASCU 



0(0,0,0) 


Ai(9,9,0) 


Bi(9,0,9) 


Oi(0,9,9) 


0(0,0,0) 


A 2 (-9 } 9,0) 


B 2 (-4,5,-ll) 


2 (3,12,-3) 


0(0,0,0) 


4,(12,3,-3) 


B 3 (7, -8,-7) 


O s (3,3,-12) 



Table 1. Coordinates of the three regular tetrahedra in Figure ^ 




Figure 1: Regular tetrahedra of side lengths 9\/2 

Along these lines and in retrospect to our work here, the following related result of Schoenberg, 
[3], who proved that a regular ra-simplex exists in Z n in the following cases and no others: 

(i) n is even and n + 1 is a square; 

(ii) n = 3 {mod 4); 

(iii) n = 1 (mod 4) and n + 1 is the sum of two squares, 

seems to open a lot more questions such as: what will be the similar corresponding parameteriza- 
tions in all these cases when regular n-simplexes exist? Also, one may ask when a regular n-simplex 
exists in Z m . 

Equilateral triangles with vertices of integer coordinates in the three dimensional space are 
numerous as one could imagine from the situation just described. One less obvious example is the 
triangle CDO with C(31, 19, 76) and D(44, 71, 11) having side lengths equal to 13\/42. Generating 
all such triangles is a natural question and we may start with one such triangle and applying the 
group of affine transformations T(x) = aO(x) + y where a £ Q, O is an orthogonal matrix with 
rational coefficients and y a vector in Q 3 . Such a transformation multiplies the side lengths with 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


£T(n) 


8 


80 


368 


1264 


3448 


7792 


16176 


30696 


54216 


90104 



Table 2. Sequence A 102698 



the factor a and so for instance the triangle OAB with side lengths y/2 cannot be transformed into 
the CDO. 

In fact, in the next section we show that the side lengths which appear from such equilateral 
triangles are of the form y/2N(<^) where £ is an Eisenstein-Jacobi integer and N(Q is its norm. The 
Eisenstein-Jacobi integers are defined as Z[o>] where u is a primitive root of unity, i.e. the complex 
numbers of the form £ = m + nu> with m, n G Z with their norm is given by N(() = m 2 + mn + n 2 . 

What makes the existence of such triangles work in space that doesn't work in two dimensions? 
We show in Proposition 12.11 that the plane containing such a triangle must have a normal vector 
n=ai+bj+ck where a, b, c are integers that satisfy the Diophantine equation 

(1) a 2 + b 2 + c 2 = 3d 2 , d<EZ. 

This equation has no non-trivial solutions, if c = for instance, according to the Gauss' charac- 
terization for the numbers that can be written as sums of two perfect squares. 

Our study of the existence of such triangles started with an American Mathematics Competition 
problem in the beginning of 2005. The problem was: 

Problem 1.1. Determine the number of equilateral triangles whose vertices have coordinates 
in the set {0, 1, 2}. 

It turns out that the answer to this question is 80. Let us introduce the notation £T(n) for the 
number of equilateral triangles whose vertices have coordinates in the set {0, 1, 2, n} for n E N. 
Some of the values of £T(n) are tabulated next: 

This sequence has been entered in the on-line Encyclopedia of Integers Sequences by Joshua 
Zucker on February 4th, 2005. The first 34 terms in this sequence were calculated by Hugo Pfoertner 
using a program in Fortran. Our hope is that the results obtained here may be used in designing 
a program that could calculate £T(n) for significantly more values of n. 

2. Planes containing equilateral triangles in 1? 

Let us denote the side lengths of an equilateral triangle AOPQ by /. We are going to discard 
translations, so we may assume that one of the vertices of such a triangle is O(0, 0, 0). If the other 
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two points, P and Q, have coordinates (x, y, z) and (u, v, w) 
before the area of AOPQ is given by 



respectively, then as we have seen 



(2) Area = i-^? = ^|OPx OQ\ = 

This implies the following simple proposition but essential in our discussion: 

Proposition 2.1. Assume the triangle AOPQ is equilateral and its vertices have integer co- 
ordinates with O the origin and I = \OP\. Then the points P and Q are contained in a plane of 
equation aa + bf3 + cry = 0, where a, b, c, and d are integers which satisfy and I 2 = 2d. 

Proof. Assume the coordinates of P and Q are denoted as before. Let us observe that 
(x — u) 2 + (y — v) 2 + [z — w) 2 = 2l 2 — 2(xu + yv + zw) which implies xu + yv + zw = £ = d G 
Z. Then using the fact that OP and OQ are contained in the plane orthogonal on the vector 
OP x OQ = ai +bj + ck with a = yw — vz, b = zu — xw and c = xv — yu the statement follows 
from ©. ■ 
The equation Q has infinitely many integer solutions besides the obvious ones a = ±d, b = ±d, 
c = ±d. For instance we can take a = —19, b = 11, c = 5 and d = 13 and the triangle OCD given 
in the Introduction has C and D in the plane {(a, (3, 7) G M 3 | — 19a + 11/3 + 57 = 0}. 

Definition 2.2. Let us introduce the setV of all these planes, i.e. {(a, (3, 7) G M 3 | aa + b(3 + 
C7 = e], such that a 2 + b 2 + c 2 = 3d 2 for some a, b,c,d,e^7, and gcd(a, b,c) = 1. 

So, if one starts with a plane n in V, picks two points of integer coordinates that belong to ir, 
which amounts to solving a simple linear Diophantine equation, the natural question is weather or 
not there exists a third point of integer coordinates, contained in tt, that completes the picture to 
an equilateral triangle (see Figure I^J). In order for the third point to exist one needs to take the first 
two points in a certain way. But if one requires only that the new point have rational coordinates 
it turns out that this is always possible and the next theorem gives a way to find the coordinates 
of the third point in terms of the given data. 

Theorem 2.3. Assume that P(u,v,w) (u,v,w G is an arbitrary point of a plane tt G V of 
normal vector (a, b, c) and passing through the origin O. Then the coordinates of a point Q(x, y, z) 



i j k 
x y z 
u v w 
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situated in it and such that the triangle AOPQ is equilateral are all rational numbers given by: 



(3) 





u 




cv — bw 


X = 




± 


2d 




2 






V 




aw — cu 


y = 




± 




2 

w 


, 2d 

bu — av 


z = 




± 


2d 


~2 






Figure 2: Plane of normal (a,b,c) 



Proof. From the geometric interpretation of the problem of finding a point Q G vr such that 
AOQP becomes equilateral one can see that there are only two points that satisfy this condition 
and just having real coordinates. We want to show that one point is given by taking the plus 
sign in all equalities in (J2J) and the other point corresponds to the minus sign in ©. We are 
going to set n = — -=(a,b,c) which is one of the two unit vectors normal to the plane it and let 

r = OP = (u,v,w). Then the cross product r x n is given by 

-» -* 1 

r x n = — i=(cv — bw, aw — cu, bu — av). 

dV3 ! 

So we observe that the solution (x, y, z) is, in fact, if written in vector notation, OQ± = \ r ± 
r x n. It is easy now to check that r and OQ± have the same norm and make a 60° angle in 
between. ■ 

3. Solutions of the Diophantine equation o? + b 2 + c 2 = 3d 2 

We have enough evidence to believe that for each (a, b, c) satisfying (QJ for some d G Z, there are 
infinitely many equilateral triangles in Z 3 that belong to a plane of normal vector (a, b, c) and the 
purpose of this paper is to determine a way to generate all these triangles. But how many planes 
do we have in VI Let us observe that if d is even, then not all of a, b, c can be odd integers, so at 
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least one of them must be even. Then the sum of the other two is a number divisible by 4 which is 
possible only if they are also even. Therefore we may reduce all numbers by a factor of two in this 
case. So, if we assume without loss of generality that gcd(a,b,c) = 1, then such solutions of ((T]) 
must have d and odd integer and then this forces that o, b and c must be all odd integers too. If 
in addition, disregarding the signs, we have a, b, c € N and a < b < c, then such a solution will be 
referred to as a primitive solution of (JJ). One can find lots of solutions of (JJ) in the following way. 



Proposition 3.1. (i) The following formulae give a three integer parameter solution of ^j: 
(4) 



a = —x\ + x\ + x\ — 2x\X2 — 2x1X3 
b = x? — x\ + x\ — 2x x x\ — 2x2X3 

2,2 2 o o , X!,X2,X 3 EZ. 

C = X\ + X2 — X3 — ZX3X1 — 2x3X2 

d = x\ + x\ + x 2 

(ii) Every primitive solution of £7J) is of the form 0]j wni/t x±,X2,X3 E Q\/fc wi/i some A; 6 N. 

Proof, (i) Let us observe that a = d — 2x\s, b = d — 2x2s, X3 = d — 2x3s where s = x± + X2 + X3. 
Then a simple calculation shows that a 2 + b 2 + c 2 = 3d 2 . In fact, for every x±,X2,X3 S Q for which 
a, b, c, d given by (@J are integers we get a solution of Q. 

If one has a solution (a, 6, c, d) of Q with d / 0, then we can introduce x\ = (d—a)t, X2 = (d—b)t, 

X3 = (d — c)t where t£l such that x\ + x 2 , + x\ = d. This gives 

9 d d lit 

t 2 - 



(d-a) 2 + (d-6) 2 + (d-c) 2 6d 2 - 2d(a + b + c) 2(3d - a - b - c) 2xi + x 2 + x 3 ' 
and so t = with s = X1+X2+X3 = \J (3d — a — b — c)/2. Every primitive solution of Q must have 
a, 6, c, and d all odd numbers as we have observed. This make k = (3d—a — b—c)/2 an integer. Since 
3d > a + b+c is equivalent to 3(a 2 + b 2 + c 2 ) - (a + b+c) 2 = (a-b) 2 + (b- c) 2 + (a- c) 2 > 0, it follows 
that fc G N. Therefore in general every solution is of the form in @ but with xi,X2,X3 G Qv^. ■ 

Example: Suppose a = 5, b = 11, c = 19 and <i = 13 then xi = 2\/2, X2 = ^a/2 and X2 = — \\[2. 

Although formulae (@J) may provide infinitely many solutions of (^Q), the following proposition is 
giving more information about its integer solutions. 

Proposition 3.2. The equation 0) has non-trivial solutions for every odd integer d > 3. 

Proof. If d = 2p + 1 for some p e Z, p > 1, then 3d 2 = 3[4p(p + 1) + 1] = 8/ + 3 which shows 
that 3q 2 = 3 (mod 8). From Gauss's Theorem about the number of representations of a number 
as a sum of three squares (Theorem 2, pp. 51 in 0), we see that the number of representations of 
3d 2 as a sum of three squares is at least 24. Here, the change of signs is counted so each solution 
actually generates eight solutions by just changing the signs. Also, the six permutations of a, b and 
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d 


(a,b,c) 


d 


(a,b,c ) 


1 


{(1,1,1)} 


9 


{(1,11,11), (5,7,13)} 


3 


{(1,1,5)} 


11 


{(1,1, 19), (5, 7, 17), (5, 13, 13)} 


5 


{(1,5,7)} 


13 


{(5, 11, 19), (7, 13, 17)} 


7 


{(1,5,11)} 


15 


{(1,7, 25), (5, 11, 25), (5, 17, 19)} 



Table 3. Primitive solutions of Q 



c get into the counting process. So, there must be at least one solution which is nontrivial since 
the solution 3d 2 = d 2 + d 2 + d 2 generate only eight solutions by using all possible change of signs 
and all permutations. ■ 
Another way of generating an infinite family of solutions of Q is to reduce it to two separate 
equations, say: 146 = 3d 2 — c 2 and a 2 + b 2 = 146. The second equation admits as solution, for 
instance, a = 11 and 6 = 5. The first equation has a particular solution d = 7 and c = 1. Then one 
can use the recurrence formulae to obtain infinitely many solutions of 146 = 3d 2 — c 2 : 

dn+i = 2d n + c n , c n +i = 3d n + 2c n for n £ N 

and d\ = 7, c\ = 1. A simple calculation shows that 3d 2 1+l — c^ +1 = 3<i„ — c 2 n so, by induction, 

{dm C-n 

) is a solution of the equation 3d 2 - c 2 = 146 for all n G N. It is easy to see that q n and c n 
are increasing sequences so this procedure generates infinitely many solutions of (^Q). 
Some of the primitive solutions of Q are included in the table below. 

As a curiosity the number of representations of d = 1003 (which will give the number of "primi- 
tive" representations of a 2 + b 2 + c 2 = 3 x 2006 2 ) is 182. 

4. Generating the beginning of the parametrization family 

The simplest solution of (^Q) is a = b = c = d = 1. We are going to introduce some more notation 
here before we give the parametrization for this case. 

Definition 4.1. For every (a,b,c), a primitive solution of (0), denote by T a b,c the set of all 
equilateral triangles with integer coordinates having the origin as one of the vertices and the other 
two lie in the plane {(a, (3, 7) G K 3 1 aa + b(3 + cy = 0} . 

It is clear now that in light of Proposition l2~l1 every equilateral triangle having integer coordinates 
after a translation, interchange of coordinates, or maybe a change of signs of some of the coordinates, 
belongs to one of the classes T^ a ,b,c- Let us introduce also the notation T for all the equilateral 
triangles in Z 3 . If we have different values for a, b and c, how many different planes can one obtain 
by permuting a, b and c in between and changing their signs? That will be 6 permutations and 
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essentially 4 change of signs (note that ax + by+cz = is the same plane as (— a)x + {— b)y+(— c)z = 
0) which gives a total of 24 such transformations. We are going to denote the group of symmetries 
of the space determined by these transformations and leave the origin fixed, by S cu b e (it is actually 
the group of symmetries of the cube). Hence we have 

(5) T= |J s(T aAc )+v. 

S 6 S cube , s(0) = O, 

a 2 + b 2 + c 2 = 3d 2 
< a < b < c, gcd(a, b, c) = 1 
a, b, c,d G Z,v G Z 3 

Theorem 4.2. Every triangle OAB G Ti,x,\ is of the form {A, B, O} = {(m, —n, n — m), (m — 
n, —m,n), (0,0,0)} for some m,n £ Z. The side lengths of AOAB m ^ n are given by 

I = \j 2{m 2 — mn + n 2 ). 

Proof. Let us assume A has coordinates (u,v,w) with u + v + w = and B(x,y,z) with 
x + y + z = 0. From © we get that x = | + y = | + z = § + ^ if we choose the 
plus signs. This choice is without loss of generality since we can interchange the roles of A and B 
if necessary. This implies x = —w, y = —u and z = —v. So, if we denote u = m, v = —n then 
w = n — m and so x = m + n, y = —m, z = —n. ■ 

We are introducing the notation N(m,n) = 2(m 2 — mn + n 2 ) for m, n G Z. For the next cases, 
d G {3,5,7}, as we have recorded in the Table |3J 3d 2 has also a unique primitive representation. 
One can use basically the same technique as in the proof of Theorem l4.2l to derive the corresponding 
parameterizations for the vertices in T a \, c (d G {3,5,7}) and the corresponding side lengths but 
for each individual set of formulae, that are given below, we had something specific to work out in 
order to get rid of denominators that naturally arise from applying @: 



d = 3, 1 = 3^N{m,n) 
Tx i 5 = {[O, (4m — 3n, m + 3n, —m), (3m + n, —3m + 4n, —n)] : m,n G Z, Z ^ 0}, 



d = 5, 1 = 5 v / iV(m,n) 
^1,5,7 = {[O, (7m — 4n, 5n, — m — 3n), (3m — 7n, 5m, —4m + n)] : m, n G Z, Z / 0}, 

c2 = 7, Z = 7y/N(m,n) 

7l,5,11 = {[O, (8m — 9n, 5m + 4n, —3m — n), (— m — 8n, 9m — 5n, —4m + 3n)] : m, n G Z, Z ^ 0}. 
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Remark: Every triangle in one particular family, s(T a i, c ) +v, is different of all the other triangles 
in other families since they live in different planes. So if we take in (J^J) only the symmetries s G S cu b e , 
s(O) = O, that give different normal vectors (two of the numbers a, b, c may be equal) then © is 
a partition of T. 

The case d = 9 is the first in which there are two essentially different primitive representations 
of 3d 2 : 3 x 9 2 = l 2 + ll 2 + ll 2 + 1 = 5 2 + 7 2 + 13 2 . The corresponding parameterizations are 
included here for completion: 



d = 9, 1 = 9y/N(m,n) 

ll = {[O, (11m — lln, 4m + 5n, —5m — An), (— lln, 9m — 4n, —9m + 5n)] :m,n£Z,l^ 0}, 
T§ : 7 t i3 = {[O, (7m + 5n, 8m — lln, —7m + 4n), (12m — 7n, —3m — 8n, —3m + 7n)] :m,n£Z,l/0}. 

To give and idea of how we obtained these parametrizations we will include the proof of the case 
d = 9, a = 5, b = 7, c = 13, that gave us the last of the above formulae. 

Proof. Assume that one of the points, P, has coordinates (m, v,w). If one solves the Diophan- 
tine equation 5u + 7v + 13w = finds that a general solution may be written as 

J w = 5u + 7t, 

\v = — lOn — 13t and t,w G Z. 
Using Theorem 12.31 we see that the coordinates of a point Q, say (x,y,z), such that AOPQ G 
75^13 must be given by ©• Switching P with Q, if necessary, we may assume that the signs in 
© could be taken all plus signs. This gives 

( _ 26u 109t 
X ~ 3 9~ ' 

13u Alt 
V ~ 3 9~' 

_ 17u 64t 

Since x = —9u — 12t — must be an integer we need to have t = 3u + 9g for some g G Z. 
Substituting we find that all other coordinates are integers: x = — 45u — 109g, y = —18m — 41g, 
z = 27n+64o. Calculating I 2 = v?+v 2 +w 2 we get I 2 = 3078m 2 + 14742u#+17658# 2 = 2x9 2 (19m 2 + 
91u#+109g 2 ). Or I 2 = 2x9 2 (19m 2 +91mo+109# 2 ) = 2x9 2 [(2u+55) 2 + (2m+55)(3m+7c/) + (3m+7c/) 2 ] 
which suggests that we could change the variables 2u + 5g = —m and 3u + 7g = n to obtain the 
statement from the theorem. If we solve this system it turns out that the solution preserves integers 
values since u = 7m + 5n and g = —3m — 2n. ■ 
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A natural question that we may ask at this point is whether or not every T a bc admits such a 
parametrization. In the next section we prove that this is indeed the case under the assumption 
that min{gcd(a, d), gcd(b, d), gcd(c,d)} = 1. We checked all the odd integers d between 1 and 401 
and this condition is satisfied for all such d's. 

5. Characterization of side lengths 

We will begin with two preliminary results. The first we just need to recall it since it is a known 
fact that can be found in number theory books mostly as an exercise or as an implicit corollary of 
more general theorems about quadratic forms or Euler's 6k + 1 theorem (see jl], pp. 568 and pQ, 
pp. 56). 

Proposition 5.1. An integer t can be written as m 2 — mn + n 2 for some m,n £ Z if and only 
if in the prime factorization of t, 2 and the primes of the form 6k — 1 appear to an even exponent. 

The next lemma is probably also known in algebraic number theory but we do not have straight 
reference for it so we are going to include a proof of it. 

Lemma 5.2. An integer t which can be written as t = 3x 2 — y 2 with x, y £ Z is the sum of two 
squares if and only if t is of the form t = 2(m 2 — mn + n 2 ) for some integers m and n. 

Proof. For necessity, by Proposition 15.11 we have to show that t is even and t/2 does not 
contain in its prime factor decomposition any of the primes 2 or those of the form 6k — 1 except to an 
even power. First, let us show that t must be even and the exponent of 2 in its prime factorization 
is odd. Since t = 3x 2 — y 2 = a 2 + b 2 implies 3a; 2 = a 2 + b 2 + y 2 we have observed that either all x, 
y, a, and b are even or all odd. 

If x, y, a and b are all even we can factor out a 2 from all these numbers and reduce the problem 
to t/4 instead of t. Applying this arguments several times one can see that t = 2 2l+1 t' with t' odd 
and I G Z. Without loss of generality we may assume that / = 0. In this case t contains only one 
power of 2 in its prime decomposition and so x, y, a and b must be all odd. 

Let us then suppose that t/2 is divisible by a prime p = 6k — 1 for some k E N. We need to 
show that the exponent of p in the prime factorization of t is even. Since p divides t = 3x 2 — y 2 
we get that 3x 2 = y 2 (mod p). If p divides x, then p divides y and so p 2 divides t which reduces 
the problem to t/p 2 . Applying this argument several times we arrive to a point when p does not 
divide x/p 1 . So, discarding an even number of p's from t, we may assume that i = 0. This implies 
that x has an inverse modulo p and then z 2 = 3 (mod p), where z = x _1 y. Using the Legendre 
symbol this says that (-) = 1. By the Law of Quadratic Reciprocity (Theorem 11.7, in [4 1) we 
see that (|) = (— 1) ( 2 ~ ~ = (— l) 3 ^" 1 . But the equation in w, w 2 = p (mod 3), is equivalent to 
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w 2 = — 1 (mod 3) which obviously has no solution in w. This implies (|) = —1 and so k has to 
be even. Therefore p = 12k! — 1 = 4j + 3 for some j £ Z. But by hypothesis, t is a sum of two 
squares and so, from Euler's characterization of those numbers, p must have an even exponent in 
the prime decomposition of t. 

For sufficiency, let us assume that t = 3x 2 — y 2 = 2(m 2 — mn + n 2 ) for some x,y,m,n G Z. 
Using Euler's characterization we have to show that if p = Ak + 3 is a prime dividing t then the 
exponent in its prime decomposition is even. If p = 3, then 3 divides y which implies x 2 — 3y' 2 = 
2{m' 2 — m'n' + n' 2 ). This is true because of Proposition l5?T1 which one has to use in both directions. 

If m' 2 —m'n' +n' 2 is not divisible by 3 then m' 2 +m'n' +n' 2 = 1 (mod 3) since all the prime factors 
of the form 6k — 1 and 2 appear to even exponents. This implies x 2 — 3y' 2 = 2(m' 2 — m'n' + n' 2 ) 
(mod 3) or x 2 = 2 (mod 3) which is a contradiction. So, ml 2 — m'n' + n' 2 must contain another 
factor of 3 and so the problem could be then reduced to t/9 instead of t. Hence 3 must have an 
even exponent in the prime decomposition of t. 

Let us assume that k = 3 j — 1 with j £ Z. Then p = 12j — 1 = 6(2j) — 1 and so these primes 
must appear to an even power in the decomposition of m 2 — mn + n 2 . The case k = 3j : + 1 
(p = 12k + 7) is not possible because that will contradict the Law of Quadratic Reciprocity: 
(I)(f) = l/(-l) 1 ^^. 



Theorem 5.3. An equilateral triangle of side lengths I and having integer coordinates in R 3 
exists, if and only if I = y / 2(m 2 — mn + n 2 ) for some integers m and n ( not both zero ) . 

Proof. The sufficiency part of the theorem is given by the triangles in T\\ \ (Theorem I4.2|) . 
For necessity let us start with an arbitrary equilateral triangle having integer coordinates and 
non zero side lengths /. Without loss of generality we may assume that one of its vertices is the 
origin. Denote the triangle as before AOPQ, with P(u,v,w) and Q(x,y,z). As we have shown 
in Proposition 12.11 we know that au + bv + wc = and ax + by + cz = for some a, b, c satisfying 
a 2 + b 2 + c 2 = 3d 2 and gcd(a, b, c) = 1. We noticed too that all a, 6, c have to be odd integers and 
so, in particular, they are all non-zero numbers. 

Then 

bv + cw \ 2 2 2 ( ffl2 + b 2 )v 2 + 2bcvw + (a 2 + c 2 )w 2 



I 2 = u 2 + v 2 + w 2 = [ ' I + v z + w 

\ a J a* 

Completing the square we have 

bcw \ 2 ,„ l0 ,a x n b 2 c 2 w 2 



aV = (a 2 + 6 2 )u 2 + 2bcvw + (a 2 + c 2 )w 2 = (3d 2 - c 2 ) L + ^jTT^ ) + ^ ~ ^ 



3d 2 - c 



2 
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or 



a 2 (3d 2 - c 2 )l 2 = [(3d 2 - c> + bcw] 2 + 3d 2 a 2 w 2 . 

This calculation shows that a 2 (3d 2 — c 2 )l 2 = m' 2 — m'n' + n' 2 where m! = (3d 2 — c 2 )v + bcw + daw 
and n' = 2daw. By Lemma 15.21 we can write 3d 2 — c 2 = a 2 + b 2 = 2(m" 2 — m"n" + n" 2 ) for some 
m",n" G Z. Hence 



I 2 = 2- 



1 m' 2 — m'n' + n' 2 



' (2a) 2 m" 2 - m"n" + n" 2 ' 
Because I 2 G Z and Proposition 15.11 we see that £ 2 = 2(m 2 — mn + n 2 ). ■ 

We include here a similar result and the last of this section which is only based on Proposition l2.ll 

Proposition 5.4. A regular tetrahedra of side lengths I and having integer coordinates in R 3 
exists, if and only if I = my2 for some m G N. 

Proof. For sufficiency, we can take the tetrahedra OPQR with P(m,0,m), Q(m,m,0) and 
R(0, m, m). 

For necessity, without loss of generality we assume the tetrahedra OPQR is regular and has all 
its coordinates integers. As before, we assume P(u,v,w) and Q(x,y,z). 




Figure 3: Regular tetrahedra 
Let E be the center of the face AOPQ. Then from Proposition 12.11 we know that = 

\ER\ 

{aM = ^(a, b, c) for some a, b,c,d G Z, I 2 = 2d. The coordinates of E are (S±£ ^). 



vWF^ ~ v 7 ^ ' ' 7 ' ' ' ' iu C Luuium« ra uiJ,ai Cl 3,3,3 

From the Pythagorean theorem one can find easily that RE = l\J\- Since OR = OE + ER, the 
coordinates of R must be given by 
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U + X 



3 V3d 



y + v t 

a, — - — ± l\ 



3-v/3d 



2 1 



3 y/M 



or 



u + x 2\/2 y + v 2\/2, 2 + w 

i „, a, — - — ± — o, — - — ± „, c 



3 3/ ' 3 3/ ' 3 3Z 

Since these coordinates are assumed to be integers we see that / = my/2. ■ 
A natural question that one may ask at this point is weather or not every family T a ^,c contains 
triangles which are faces of regular tetrahedra with integer coordinates. We believe that every 
triangle in T a ^ fi is the face of such a tetrahedra as long as its sides, in light of Theorem 15.31 are 



of the form I = \/2N(m, n) with N(m, n) a perfect square. We leave this conjecture for further 
study. To find values of m, n £ Z such that N(m, n) 



m — ran + n is a perfect square, of course 
one can accomplish this in the trivial way, by taking m = or n = but there are also infinitely 
many non-trivial solutions as one can see from Proposition 15.11 

6. A MORE GENERAL PARAMETRIZATION 

Our construction depends on a particular solution, (r, s) G Z 2 , of the equation: 

(6) 2(a 2 + b 2 ) =s 2 + 3r 2 . 

As before let us assume that a, b, c and q are integers satisfying a 2 + b 2 + c 2 = 3d 2 with d an odd 
positive integer and gcd(a, b, c) = 1. By Lemma l5~2*l we see that 3d 2 — c 2 = a 2 + b 2 = 2(f 2 — fg + g 2 ) 
for some f,g e Z and so 2(a 2 + b 2 ) = (2f - g) 2 + 3g 2 which says that equation © has always an 
integer solution. 

Theorem 6.1. Let a, b, c, d be odd positive integers such that a 2 + b 2 + c 2 = 3d 2 , a < b < c 
and gcd(d, c) = 1. Then T a ,b,c = { AOPQ\ m,n G Z} where the points P(u, v, w) and Q(x, y, z) are 
given by 



(7) 
with 



u = m u m — n u n, 
v = m v m — n v n, 
w = m w m — n w n, 



and 



m x m — n x n, 



z = m z m — n z n, 



n x = —(rac + dbs)/q 
n y = (das — bcr) / q 



n. 



(8) 



m x = — ^[db(3r + s) + acir — s)]/q, 
m y = ^[da(3r + s) — bc(r — s)]/q, 
m z = (r - s)/2, 
and 

m u = —(rac + dbs)/q, n u = — ^[db(s — 3r) + ac(r + s)]/q 
m v = (das — rbc)/q, n v = ^[da(s — 3r) — bc(r + s)]/q 
m w = r, n w = (r + s)/2 
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= a 2 + b 2 and (r, s) is a suitable solution of JZS)). 



In order to show Proposition 16. II we need the following lemma. 

Lemma 6.2. Suppose A and B are integers in such a way A 2 + 3B 2 is divisible by q where 2q 
can be written as s' 2 + 3r' 2 with r' ,s' G Z. Then there exist a writing 2q = s 2 + 3r 2 , r, s 6 Z, such 
that 

Ar + Bs = (mod2q), 

and 

As - 3Br = (mod 2q). 

Proof. Let us observe that every number of the form A 2 + 3B 2 is an Eisenstein- Jacobi integer 
since A 2 + 3B 2 = (A + B) 2 - (A + B)(2B) + (2B) 2 . Conversely if n is even then m 2 - ran + n 2 = 
(m — n/2) 2 + 3{n/2) 2 and since m 2 — ran + n 2 = [n — ra) 2 — (n — ra)n + n 2 we see that every 
Eisenstein-Jacobi integer is of the form A 2 + 3B 2 for some A, B S Z. Using Proposition 15. II we can 
write 

A 2 + 3B 2 = 2 2a (l[p t ) 2 l[p j , a £ N, 

teT jeJ 

and 



2q = 2 2 ?{\[p t ) 2 Y[ Pj , 1 < & < a, T' C T, J' C J, 

te/' is J' 

with pt prime of the form 6k — 1 for i 6 T and pj prime of the form 6k + 1 or equal to 3 for all 
j G J. From Euler's 6k + I theorem, for each j G J we can write = (m-j + nj^/3i)(m,j — nj 
and we make the choice of mj and nj in Z in such a way yl + B^3i = s (mj + n^i) using the 

prime factorization in T\\[3i\ of A + B\[3i and /i = 2° F] p<. 

Then we take r and s such that s + r\[3i = h' F] (m^ — nj\/3i) with /i' = 2^ F] p^. Notice that 

jeJ' t&T' 

2q = (s + rV3i)((s - rV3i) = s 2 + 3r 2 and (A + By/3i)(s + rV3i) = 2q(A' + 5'\/3i). Identifying 
the coefficients we get As + 3Br = 2qA' and ^4r + Bs = 2qB' and the conclusion of our lemma 
follows from this. ■ 
To return to the proof of Theorem 16. II we begin with the next proposition. 

Proposition 6.3. For some particular solution (r,s) £ Z 2 of @) all of the coefficients m u , 
ra v j ra w j nu, ny, n w , ra x , ray, raz, rix,ny, n z %n JSp are integers. 
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Proof. One can check that 
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0) 



m 2 + m 2 +m 2 z = n 2 + n 2 + n 2 = m\ + m 2 + m\ = n 2 u + n 2 +n 2 w = 2d 2 , 
m x n x + m y n y + m z n z = m u n u + m v n v + m w n w = d 2 , 

am x + bm y + cm z = am u + bm v + cm w = an x + bn y + cn z = an u + bn v + cn w = 



These identities insures that the points P(u, v, z) and Q(x, y, z) are in the plane of normal vector 
(a, b, c) and containing the origin, the AOPQ is equilateral for every values of m, n and its side 
lengths are I = dy / 2(m 2 — mn + n 2 ). These calculations are tedious and so we are not going to 
include them here. The only ingredients that are used in establishing all these identities are the 
two relations between a, b, c, d, r and s. 

From © we see that r and s have to be of the same parity. Then, it is clear that m z , n z , m w , 
and n w are all integers. Because the equalities in © are satisfied it suffices to show that m x , n x , 
m u , and n u are integers for some choice of (r, s) solution of ©• To show that n x is an integer we 
need to show that q divides N = rac + dbs. 

Let us observe that c 2 = 3d 2 and a 2 = —b 2 (mod q). Multiplying together these two congruencies 
we obtain (ac) 2 + 3(db) 2 = (mod q). Using Lemma 16.21 we see that N is divisible by 2q for some 
choice of r and s as in ©. So, n x £ Z. 

Next we want to show that m x is an integer. First let us observe that if M = 3dbr — acs we 
can apply the second part of Lemma 16.21 to conclude that 2q divides M also. Hence 2q divides 
M + N = db(3r + s) + ac(r — s) and so, m x is an integer. Because m x + n u = n x and m u = n x it 
follows that n u and m u are also integers. ■ 

Remark: Let us observe that the parametric formulae (J7| and (JSJ) exist under no extra assump- 
tion on a, b and c. The questions is if every triangle in T a ,b,c ls given by these formulae. 

So, with these preparations we can return to prove this under Theorem 16. li s assumption. 

Proof. We start with a triangle in T a ^ c say AOPQ with the notation as before. We know 
that P(uo,vo,wo) and Q(xo, yo, Zo) belong to the plane of equation aa + b(3 + C7 = and by 
Theorem 12.31 we see that the coordinates of P and Q should satisfy ©• Hence, using the same 
notation, cvq — bwo, clwq — cuq and buQ — awo are divisible by d. A relatively simple calculation 
shows that m v n w — m w n v = ad, m w n u — m u n w = bd and m u n v — m v n u = cd. We would like to 
solve the following system in m and n: 



(10) 



= rn u m — n u n, 
vo = m v m — n v n, 
wq = m w m — n w n. 
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The equalities © and the fact that (uq,vq, wq) is in the plane aa + b(3 + cy = insures that (fTUl) 

has a unique real solution in m and n. We want to show that this solution is in fact an integer 

solution. The value of n can be solved from each pair of these equations to get 

_ vom w - wom v _ w m u - u m w _ u Q m v - v m u 
ad bd cd 

Since gcd(a,b,c) = 1 we can find integers a', b' , d such that aa' + W + cd = 1. Hence the above 

sequence of equalities gives 

_ a'(vom w - w Q m v ) + b'(wom u - u m w ) + c'(u m v - v m u ) 
H ~ d 
So, it suffices to show that d divides vom w — wom v , wom u — uom w and uom v — VQm u . 

Next, let us calculate for instance vom w — wom v in more detail: 

das — rbc v$qr — (das — rbc)wQ 

v m w - w m v = v r w = = 

q q 

—daswo + vo(3d 2 — c 2 )r + rbcwo c(bwo — vqc)v + 3rvod 2 — daswo 

q q 

From Theorem 12.11 we see that bwo — vqc = ±d(2xo — uq)- Hence, 

d[±c(2xo — uo)r + 3rvod — aswo] 
v m w - w m v = . 

q 

This shows that d divides vom w — wom v provided that gcd(q, d) = 1. This is true since gcd(d, c) = 
1 implies gcd(q, d) = 1. Similar calculations show that d divides wom u — uom w and uom v — VQm u . 
Hence n must be an integer. Similarly one shows that m is an integer. By replacing P with Q if 
necessary all the equalities in (JJJ) have to hold true by Theorem 12. II ■ 

It is natural to ask if the Diophantine equation a 2 + b 2 + c 2 = 3d 2 has any solution which satisfies 
gcd(a,b,c) = 1, gcd(a,d) > 1, gcd(b,d) > 1 and gcd(c,d) > 1. Certainly the parametrizations (JJJ) 
and Q are the most general that one can hope for in case the answer to the above question is 
negative. Instead of settling this question maybe one can say how rare a solution like that can be 
and find in the process an asymptotic formula for £T(n). Let us use a similar notation here for 
the number of regular tetrahedra with coordinates in the set {0, 1,2, ...,n}: HT(n). Is there any 
simple relation between £T(n) and lZT(n)7 What is an asymptotic formula for lZT(n)l 
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